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Abstract 

The range corrections to the universal properties and structure of two-neutron halo nuclei are 
investigated within an effective quantum mechanics framework. Treating the nucleus as an effective 
three-body system, we make a systematic improvement upon previous calculations by calculating 
the linear range corrections at next-to- leading order. Since the effective ranges for the neutron- 
core interactions are not known, we estimate the effective range to be set by the inverse of the 
pion mass. We investigate the possibility of excited Efimov states in two- neutron halo nuclei 
and calculate their mean square radii to next-to-leading order. We find that the effective range 
corrections are generally small and the leading order predictions are very robust. 
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I. INTRODUCTION 



There is a considerable interest in the universal properties of physical systems with large 
scattering lengths The low-energy behavior of a few-body system with short-range in- 

teractions and scattering length \a\ much larger than the range R of the underlying two-body 
interaction can be described in a model-independent way. This is achieved by formulating 
an effective theory that takes advantage of the separation of scales inherent in the system, 
thereby using the ratio R/\a\ as a small expansion parameter. The theoretical uncertainty 
can then be systematically improved by including increasingly higher order terms in the 
expansion. These techniques can readily be applied to halo nuclei: a special class of nuclear 
systems consisting of a tightly bound core and a "halo" of lightly bound nucleons (5j-[8j]. 
There is a growing history of studies which exploit the separation of scales found in halo 
nuclei to describe the system in an effective theory, assuming the core to be a structureless 
particle and treating the nucleus as a few-body system of core and valence nucleons (9 16 



In previous work [16|, we explored the universal properties and structure of 2n halo 
nuclei to leading order (LO) in the expansion in R/\a\ by describing the halo nuclei as an 
effective three-body system consisting of a core and two lightly bound valence neutrons. We 
constructed an effective potential made up of two-body contact interactions, and showed 
that this approach is equivalent to an effective field theory (EFT) framework. Our main 
focus was the possibility of such three-body systems to display the universal Efimov effect 



17|| . Assuming the nuclear ground state to be an Efimov state, we confirmed earlier model 
studies showing that 20 C is the only 2n halo nucleus candidate to display an excited Efimov 
state (9], [To| if the 19 C binding energy is less than 165 keV. We also calculated the properties 
of this state. Moreover, we studied the structure of other known halo nuclei, calculating the 
matter density form factors and mean square radii (See also Ref. [14J for a previous study). 

In this paper, we improve upon our earlier calculation by including the momentum de- 
pendent next-to-leading order (NLO) term in the short-range effective potential. While the 
scattering of two particles at sufficiently low energy is determined to leading order by only 
the scattering length a, at NLO, the effective range tq of the two-body interaction comes 
into play. Including the linear range correction, we should be able to predict low-energy 
observables up to errors of order (r /a) 2 ~ (M[ ow /M high ) 2 . 

Effective field theories with contact interactions have already been used to study the NLO 
range corrections to various few-body systems with large scattering lengths, mostly focused 
on the three-nucleon system. The range corrections to S-wave neutron-deuteron scatte ring 
in the doublet channel were calculated in perturbation theory in Ref. (l8j. In Refs. (lOl 



2 If . low-energy three-body observables were calculated up to next-to- next-to-leading order 



(N2LO) using a formalism that partially resums higher-order range effects. An analysis of 
three- and four-body observables to NLO within the resonating group model was carried out 
in [22| . For a more exhaustive list of such higher-order studies in the three-nucleon system, 
see Refs. 0, 0] and references therein. 

There has also been interest lately in calculating the effective range corrections to the 



Efimov effect for three identical bosons. Th0gersen et al. [23| focused on the range correc- 
tions for ultracold atoms caught in an optical trap. Platter, Ji, and Phillips have explored 
the corrections to the Efimov spectrum linear in ro using once-subtracted momentum-space 
integral equations derived from effective field theory [24). They showed that discrete scale 
invariance relates the relative corrections for different Efimov states and the bound state 
spectrum in the unitary limit is unaffected by the linear range corrections. 
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In this work, we focus on the NLO correction to the properties of two-neutron (2n) 
halo nuclei arising from a non-zero effective range in the neutron-neutron and neutron-core 
interactions. In the next section, we discuss the theoretical framework for this extension. 
Our results are presented in Section IIII1 We discuss the impact of this correction on the 
possibility of known 2n halo nuclei to have an excited Efimov state and on their binding 
energies. Moreover, we calculate the shift in the mean square radii of various halo nuclei 
due to the non-zero effective range. Our conclusions are given in Section IIVI Finally, some 
details on the derivation of the Wigner bound and the calculation of the matter form factors 
are given in the appendices. 



II. THEORETICAL FRAMEWORK 



We want to extend the effective quantum mechanics framework used in Refs. [16J, [25| to 
next-to-leading order (NLO). This formulation is equivalent to using field-theoretic language 
for the problem at hand. The short-range interactions characteristic of halo nuclei are 
described using an effective interaction potential. Using the NLO effective potential, the 
low-energy behavior of the system will be reproduced up to errors proportional to the low- 
momentum scale Mi ow over the high-momentum scale squared. At this order, we 
require two coupling parameters C and C 2 tuned to reproduce the scattering length a and 
effective range r of the two-body interaction. The theory is valid up to a momentum scale 
Mhigh at which the errors become of order one. 

For the large separation of scales involved in halo nuclei, zero-range interactions can be 
used in constructing the effective interaction potential. This leads to a separable potential 
made up of contact interactions in a momentum expansion. The two-body S-wave NLO 
potential can be written as 

(p I Veff I v) = g(p)g(p)(c + c 2 (p 2 + P ' 2 ) + ...), (1) 

where the dots indicate higher order interactions and g(p) is a regulator function (sometimes 
called the form factor). Of course, the low-energy observables must be independent of the 
regularization scheme, and one can choose the scheme most convenient for calculations. 
For this study, we use a strong cutoff regularization, in which the effective potential is 
set to zero, V e ff = 0, for momenta p,p' > A. This corresponds to the regulator function 
g(p) = 6 (A — p) and allows for a simple inclusion of a three-body force term, as described 
below. A natural choice for the value of A is A ~ M^h, but observables are independent of 
A after renormalization, up to higher order corrections which scale with powers of 1/A. 

At this point it is worth commenting on the antisymmetry of the full wave function 
with respect to the exchange of nucleons in the core and the halo neutrons. In microscopic 
cluster models, such as the resonating group model [26] or fermionic molecular dynamics 



271 ]. the wave function is antisymmetric with respect to exchange of all nucleons. In our 
halo effective theory the wave function is only antisymmetric with respect to exchange of 
the halo nucleons since the nucleons in the core are not resolved. The error involved in this 
approximation is governed by the ratio of the binding momentum of the halo nucleons to 
the binding momentum of the core, i.e. Mi ow /Mhi g h- Within the domain of validity of the 
effective theory such effects are subsumed in the effective range parameters of the neutron- 
core interaction and can be treated in perturbation theory. This can be understood as 
follows: The exchange of a nucleon from the core and a halo nucleon can only give a sizeable 
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contribution if there is a spatial overlap between the wave function of the core and the wave 
function of the halo nucleon. The overlap of the wave functions, however, is determined by 
the ratio M\ ow / 'Mhigh- If this ratio is of the order 1/3, one can expect the antisymmetrization 
effects to be 30% at leading order and 10% at next-to-leading order. For a microscopic study 
of antisymmetrization effects in the neutron-alpha system, see Ref. (28| . 

The potential ([T]) is then used in the solution of the three-body Faddeev equations in 
terms of the spectator functions Fi(q), which represent the dynamics of the core (i = c) 
and the halo neutron (i = n). To find the bound state of a halo nucleus composed of two 
valence neutrons and a core with spin zero, the resulting coupled integral equations are 
simply a generalization of the three-boson equation 1 (see 25|, |29] and references therein). In 
the remainder of the paper, we use units such that h — c — 1. For convenience, we set the 
nucleon mass m — 1. The resulting integral equations for the spectator functions F n (q) and 
F c (q) and the bound state energy B 3 > are: 



F n {q) 



2 Jo 



'„'2 



dqq 



G n (q,q';B 3 ) + ^) t n (q' ; B 3 ) F n (q') 



+ [G c (q,q';B 3 ) + 



g(A) 
A 2 



A 2 



t c (q';B 3 )F c (q') 



(2) 



Fc(q) 



n 



G c (q',q;B 3 ) + 



H(A) \ 
A 2 J 



t n {q'-B 3 ) F n (q') 



(3) 



The functions G c and G n arise from the free three-body propagator. They are given by 



A 

G n (q,q';B 3 ) = — In 
qq' 



— l B 3 + q +q + a+iQQ 



2A 
A+ 



_£> 3 + q 2 + q a _ _L_ qq , 



(4) 



G c {q,q;B 3 ) = — In 



fi 3 + q 2 + ^-q' 2 - 99' 



(5) 



where A is the number of nucleons in the core. 

The effects of the interactions are contained in the two-body T-matrices which are ob- 
tained by solving the Lippmann-Schwinger equation for the neutron-neutron and neutron- 
core interactions with the effective potential in Eq. <§Q). In the kinematics of Eqs. (J21EJ), we 
have for the neutron-core T-matrix: 



t n (q';B 3 



A + l 

ttA 



1 - E n (q'] B 3 ) h nc (x n 



i ~ r~ 2 

— f(x nc ) + V E n (q'] B 3 ) - arctan(xr 



+h nc (x nc ) E n (q'; B 3 ) -A - J E n (q>; B 3 ) - arctan(x n ; 

7T v 7T 



-i -1 



(6) 



1 In fact, the equations are the same for any bound three-body system of two identical particles and a core 
with spin zero, which interact through the pair- wise zero-range potentials given in Eq. (TT]). 
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with 



2 A f n A + 2 



and for the neutron-neutron T-matrix: 



2 / 

t c (g'; B 3 ) = -II- £ c (g'; .63) h nn (x nn ) 

7T V 



1 2 

/(£nn) + V ^c(?'; 5 3 ) - arctan(x c ) 



with 



+^n(5„ n ) ^c(?'; -B 3 ) ( -A - a/ E c (q'; B 3 ) - arctan(a; c 

7T V 7T 



(8) 



E c (g'; J B 3 ) = J B3 + ^-g /i . (9) 

A single index n, c indicates the spectator particle and a nn , r nn , and a nc , r nc are the n-n 
and n-c scattering lengths and effective ranges, respectively. For brevity we have defined 
the variables: 

%i = — = and x n j = — j , where i = n,c. (10) 

Here, i£ n j is the two-body bound (virtual) state energy calculated from the pole of the two- 
body T-matrix at NLO in the limit A — > 00, with [i the corresponding reduced mass. In this 
limit, which implies a cutoff much larger than the momentum scales involved in the problem: 
A 3> l/\a\, yEi, V E ni , these T- matrices reproduce the usual effective range expansion at 
NLO. Therefore, the pole position can be found through the relation: 

E m =(l- y/l-2r ni /a m >) 2 /(2/xr^). (11) 

At leading order (r n j = 0), this expression reduces to the universal formula for the two- 
body binding energy for systems with a large positive scattering length: E ni = l/(2/xa^). 
Moreover, we have simplified the expressions for the two-body T-matrices (El E]) by defining 
the functions 

f(x) = - ( arctan (£) + -) (12) 

7T \ X J 

which approaches one as x becomes large and 

h (x ) - 1 2 ~ n f^ Ar ^ 2 % - nc fi 3 ) 

which vanishes for large x n i. The term proportional to h n i(x n j) A generates the effective 
range term in the two-body T-matrices (0[8]). 

Finally, the dimensionless function H(A) is a short-range three-body force. This term 
was included in the integral equations (0 [3]) in analogy to the three-boson case which our 
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equations reproduce in the limit of identical particles. It corresponds to a three-body contact 
interaction which is required for consistent renormalization For a given value of A, 

the three-body term H can be tuned to reproduce a given three-body observable and then 
other low-energy observables can be predicted using the same value H . The running of the 
three-body force H(A) with the cutoff A is governed by a limit cycle. As a consequence, 
it is always possible to find a value of the cutoff A at which the three-body force vanishes. 
Hence, at leading order A can simply be used as the three-body parameter in practical 
calculations [31]. The correct renormalization in the effective theory was explicitly verified 
25[. Apart from the three-body observable, our input parameters at NLO are the n-c and 



n-n scattering lengths (or energies) and the corresponding effective ranges. The three-body 
binding energies are then given by the values of B3 for which the coupled integral equations 
(121 E]) have a nontrivial solution. 

For most halo nuclei, the n-c scattering length is not as well known as the two-body 
bound (virtual) state energy. Therefore, we will generally use the two-body energy E nc as 
an input parameter. From E nc and r nc , we can calculate the corresponding a nc to be used 
in the two-body T-matrices through Eq. The n-c effective ranges for most halo nuclei 

are not known experimentally. However, the range of the interaction between two nucleons 
is determined by the exchange of pions. We therefore use the inverse of the pion mass, 
m w = 140 MeV, to estimate the n-c effective range, r nc l/m w = 1.4 fm. This is equivalent 
to setting the effective range to the natural low-energy length scale of the system. The value 
of the n-n effective range is known experimentally: r nn = (2.75 ± 0.11) fm [32{|. Its size 
is comparable to the inverse of the pion mass scale. In our numerical calculations we use 
r nn r nc m 1/rrin = 1.4 fm. This choice is sufficient for our purpose of estimating the size 
of the range corrections and testing the robustness of the leading order results. Moreover, 
it allows a slightly larger range of cutoffs to be consistent with the Wigner bound. 

The Wigner bound is based on general arguments and constrains the value of the effective 
range that can be generated by short-range potentials. It has been shown that the renormal- 
ization of an effective potential made up of contact interactions in a momentum expansion 
can only be performed if certain constraints are placed on the effective range [33T |. More 
specifically, the effective range must be negative in the zero-range limit of the potential. In 
our case, this corresponds to the limit A — > 00 since a finite cutoff generates a finite range. 
This constraint follows directly from a general bound on the derivatives of the phase shifts 
derived first by Wigner in 1954 [34]]. Wigner derived this bound from the fundamental prin- 
ciples of causality and unitarity, holding that a scattered wave cannot leave the scatterer 
before the incident wave has reached it. Phillips and Cohen derived the relation between the 
Wigner bound on the derivatives of phase shifts with the constraint on the effective range 
for short-range potentials [35]. This bound applies even if the potential does not vanish 
exactly outside some finite radius but merely decreases fast enough for the wave function 
to approach the asymptotic solution sufficiently quickly. The Wigner bound for zero-range 
potentials was shown to hold true no matter how many terms in the momentum expansion 
are included in the potential [33T ]. A very recent study has generalized the Wigner Bound 
to arbitrary dimension and angular momentum j36|. An explicit derivation of the Wigner 
bound constraint for our effective potential is given in Appendix |X] In our numerical results 
in the next section, the cutoff is always chosen in accordance with the Wigner bound. 
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III. EFFECTIVE RANGE CORRECTIONS TO 2n HALO NUCLEI 



In this section we present numerical results for the effective range corrections to low- 
energy observables of three-body halo nuclei composed of a core and two valence neutrons. 
We begin with the range corrections to the possibility of the Efimov effect in In halo nuclei, 
as well as to the binding energy of a possible Efimov excited state in 20 C. After that, we 
estimate and discuss effective range corrections to the mean square radii results. 



A. Range Corrections to In Halo Nuclei Binding Energies 



In Ref. [16j, we have calculated boundary curves representing the existence of an excited 
Efimov state for various values of the core mass in the parametric region defined by the 
ratios (E nc / B^) 1 ^ 2 versus (E nn / B^) 1 / 2 , where E nc and E nn are the n-c and n-n two-body 
energies related to the S-wave scattering length to LO. An analogous study was carried 
out previously in Ref. [Io| within the renormalized zero-range model and motivated our 
investigation. It should be noted, however, that this type of analysis implicitly assumes that 
the ground state of a halo nucleus is an Efimov state. More general scenarios where only 
the excited state has universal character cannot be excluded. The curve itself is built up 
of the points for which the B^ +1 ^ binding energy is equal to the scattering threshold. This 
is equivalent to finding the critical scattering lengths of the three-body system. In similar 
studies looking at the range corrections to the Efimov effect in the three-boson system 
HI S3, it was observed that including a non-zero effective range shifted the positions of 



the critical scattering lengths. However, for all values of effective ranges r nc and r nn within 
the range of validity of our effective theory {r ni <C \a n i\, i = n, c), we find these shifts to be 
extremely small such that the boundary curves calculated with a non-zero effective range 
are indistiguishable from those at LO (See Fig. 2 of Ref. (HI). 

As a consequence, the main conclusion from [16| is only slightly changed by the NLO 
range corrections: if the n- 18 C bound state energy satisfies E nc < 161 keV, 20 C will have an 
excited Efimov state. In the following, we investigate the dependence on the n- 18 C bound 
state energy in more detail. In Ref. [16j, we calculated the value of the possible excited state 
energies at LO as a function of E nc , using the standard value for a nn = (—18.7 ± 0.6) fm 



371 ]. and fixing the cutoff A to reproduce the experimental value of the ground state energy 



B^ = 3506 keV [38|, |39|. Now using the estimate r$ = r nn = r nc « l/m n = 1.4 fm for the 



effective ranges along with the standard value for a nn given above, we calculate the excited 

5 3 



state energies as a function of E nc to NLO, by fixing A and tuning the three-body term 



H to reproduce the experimental value of the ground state energy £< 0) = 3506 keV. For our 
calculations we use the cutoff A = 580 MeV. The NLO result is compared to the LO result 
in Fig. [U where the solid black line is the excited state energy to LO, the dotted-dashed red 
line to NLO, and the dashed line represents the scattering threshold. The inset graph shows 
the excited state energy relative to the scattering threshold. Whereas the excited state in 
the LO calculation exists when E nc < 165 keV, the state to NLO only exists when E nc < 161 
keV. For larger values of E nc , the 20 C system moves across the scattering threshold, and the 
excited Efimov state disappears. We also see that the effective range corrections only lead 
to a recognizable shift toward smaller binding energies for E nc < 90 keV. For larger values 
of E nc there is a region where the effective range corrections lead to a shift towards larger 
binding energies. However, this shift is always very small, no greater than « 0.5 keV over 
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200 




FIG. 1: Binding energy of the C excited Efimov state as a function of the n- C bound state 
energy to LO (solid black line) and NLO (dotted-dashed red line) with NLO error bands. The 
dashed line represents the scattering threshold which is given by = E nc . The inset shows the 
excited state energy relative to the scattering threshold. 



the whole range of possible E nc values shown, and always smaller than the NLO error bands. 

We have estimated this NLO error using the theoretical uncertainty of our effective poten- 
tial. The uncertainty in binding energies calculated using the two-body effective potential 
of Eq. ([T]) is ~ (r /a) 2 . For effective ranges much smaller than the scattering length we 
know that l/a^ c ~ 2fi nc E nc . Therefore, we estimate the theoretical uncertainty in the bind- 
ing energy of the excited state using the sum of the uncertainties from the n-n and n- ls C 



interactions: AB^/B^' « (r nn /a nn ) 2 + rl c (2fi nc E nc ). 



4" 



The TUNL nuclear data evaluation from 1995 |38|, [39J quotes the n- C bound state 
energy as E nn = (162 ± 112) keV. This value is based on experiments using strong breakup 



reactions |40l443j and would leave the possibility of an Efimov excited state in 20 C open. 
The more recent atomic mass evaluation AME2003 [441 ]. however, quotes E nc = (580 ± 100) 
keV which excludes an excited Efimov state. This evaluation includes a recent measurement 
of E nc using Coulomb dissociation of 19 C that found E nc = (530 ± 130) keV j45|. Typel 
and Baur showed that the value of Ref. 45j is also consistent the distribution of the B(E1) 
continuum strength j46|. A review of the experiments was given in (47| . 



B. Range Corrections to In Halo Nuclei Mean Square Radii 

In this subsection, we will calculate the NLO corrections to the mean square radii of 2n 
halo nuclei. We will briefly review the necessary formalism and then present our results. A 
more detailed description of the formalism is given in Refs. 0, 12911 . For the calculation of 



the one- and two-body matter density form factors and mean square radii, we need the full 
three-body wave function, which can be reconstructed from the Faddeev spectator functions 
F n and F c . The relevant formulae are given in Appendix IB1 

The mean square radii are then calculated from the matter density form factors in the 
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low momentum transfer region, where the slope of the form factor determines the mean 
square radius (r 2 ): 

F{k 2 ) = 1 - \h? (r 2 ) + . . . . (14) 

Of course, the mean square radius acquired depends on the choice of one- or two-body form 
factor. These matter density form factors are calculated from the full wave function of the 
three-body system through Eqs. (IB81 IB91 IBlOj) . The slope of the two-body form factor 



Fniik 2 ), where i — n,c, will give the mean square distance between the two particles in the 
chosen two-body subsystem, either (r 2 n ) or (r 2 c ). Analogously, the slope of the one-body 
form factor ^(k 2 ) will give the mean square distance of the spectator particle from the 
center of mass of the two-body subsystem, either (r 2 _ nn ) or (r 2 _ nc ). However, it is more 
useful to calculate the distance of the individual particles from the center of mass of the 
three-body bound state. If bi is the slope of the one-body form factor Fiik 2 ) at k 2 = 0, the 
mean square radius of one of the bodies i from the three-body center of mass is given by: 



wh^H 1 -^;- « i5) 

where is the mass of the desired constituent, % = n,c, and m n and m c are the neutron 
and core masses, respectively. 

We have extracted the radii by fitting a polynomial in k 2 to the form factor results for 
small k 2 . We have used polynomials of varying degree up to 5th order in k 2 in order to 
verify the stability and convergence of the fit. We have found a satisfactory stability in the 
slope when fitting to a polynomial to the fourth order in k 2 , up to a value of k 2 at which 
the form factor has dropped less than 10 percent. 

The extracted radii for known halo nuclei to both LO and NLO are shown in Table [B 



We show only a selection of the LO results found in Ref. [16j, to give a general overview of 



NLO corrections to the LO results. As input we have used the standard value of the n-n 



scattering length, a nn = (—18.7 ±0.6) fm 37J, along with the experimental values of the n-c 



two-body energies E nc shown in the third column of Table d (negative values correspond to 
virtual energies). The effective range of both two-body subsystems is shown in the fourth 
column, indicating which results refer to LO (ro = fm) or NLO (vq = 1.4 fm) calculations. 
As before, we use the inverse of the pion mass to estimate the effective range of both the n-n 
and the n-c interactions for all 2n halo nuclei. As a three-body input, the three-body term is 
tuned to reproduce the experimental ground state binding energy £?g shown in the second 
column of Table [I] These experimental values for the two-body and three-body energies 
are taken from the most recent results of the "Nuclear Data Evaluation Project" of TUNL 
H|, except where otherwise noted. The calculations are performed at a fixed cutoff A. Due 



to the Wigner bound, some care must be taken in choosing A. For the range of scattering 
lengths, or corresponding two-body energies, which we wish to explore, the Wigner bound 
leads to a maximum value for the cutoff for most of the 2n halo nuclei of approximately 
A < 630 MeV. For our calculations, we use a cutoff A = 580 MeV for all nuclei. This value 
is significantly larger than the typical momenta in halo nuclei which are well below the pion 
mass. 

The theoretical error for the NLO calculation is again estimated by the uncertainty of 
the two-body effective potential, Eq. (JTJ), which is of order (r /a) 2 , where r is the effective 
range of the interaction, and a the scattering length. The uncertainty in the radii is then 
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2.3±0.3 


0.19±0.02 
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3506 


60 


0.0 


2.8±0.2 


2.3±0.2 


2.2±0.2 


0.18±0.01 




3506 


60 


1.4 


2.84±0.03 


2.41±0.03 


2.28±0.03 


0.192±0.002 




65.0±6.8 


60 


0.0 


42±3 


38±3 


41±3 


2.2±0.2 


20q* 


64.9±0.7 


60 


1.4 


43.2±0.5 


38.7±0.4 


42.9±0.5 


2.26±0.02 



TABLE I: Various mean square radii of different halo nuclei. The second two columns show the 
input values for the three-body ground state energy and the two-body n-c energy (negative values 
corresponding to virtual energies), respectively, as given by 39J, except where otherwise noted. The 
fourth column shows the input value for both two-body effective ranges, related to LO (rn = 0.0 
fm) or NLO (rn = 1.4 fm) calculations. The rows marked by 20 C* show the results for the excited 
Efimov state of 20 C, with binding energy displayed in the second column, which is found above the 
ground state (B3 = 3506 keV). 



calculated exactly as was done in the previous subsection, using the sum of the uncertainties 
from the n-n and n-c interactions: A(r 2 )/(r 2 ) « {r nn /a nn ) 2 + r 2 c (2fi nc E nc ). We stress that 
our NLO results as well as their errors are based on our estimate of the effective range, 
r n = 1.4 fm. They do not include the uncertainty in the estimate of r and therefore must 
be interpreted with some care. 

We will now discuss the results for the NLO corrections to the mean square radii due 
to the effective range of the interactions. For the Borromean halo nuclei n Li and 14 Be, in 
which all two-body subsystems are unbound, the general tendency is for a positive effective 
range to shift all mean square radii to smaller values. The only exception is the case of n Li 
using the central value of the n-c energy E nc = (—25 ± 15) keV [39]. This difference is due 
to the fact that this value of the virtual energy is very close to the resonant limit, E nc = 0.0, 
while the competing value E nc = (—800 ± 250) keV 48[ is much larger in comparison. This 
can be seen more clearly if we look at the mean square radii over a range of E nc values. 
Using the central value of the three-body binding energy as input, ^ 0) = 247 keV, the NLO 
results are plotted in comparison to the LO values in Fig. [2], with error bands estimated 
from the theoretical uncertainty, as described above. The solid black lines represent the LO 
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FIG. 2: The various mean square radii for n Li as a function of the n- 9 Li two-body energy with error 
bands from the theoretical uncertainty. The LO results represented by the solid black lines, with 
error bands represented by dashed lines (as calculated in Ref. la]). The NLO results represented by 
the red dotted-dashed lines with solid error bands. As input, the n-n scattering length a nn = —18.7 
fm, the three-body binding energy = 247 keV, and for NLO results the effective range tq = 1.4 
fm were used. 



results, with error bands represented by the dashed lines Ty]. The dotted-dashed red lines 
represent the NLO results, with solid colored error bands. We see that the mean square 
radii to NLO show the same behavior with changing E nc as the LO case, with a decrease 
in the magnitude of E nc leading to an increase in size, with a more rapid increase as the 
resonant limit is approached. However, for the NLO results, the increase is more rapid when 
E nc approaches zero, such that below some value of \E nc \, a positive effective range shifts 
the mean square radii to larger values. Otherwise, a positive effective range always shifts 
the mean square radii to smaller values. One can see that this shift is more drastic as the 
magnitude of \E nc \ grows larger. This is most likely due to the fact that for these large 
values of \E nc \, the scattering length becomes smaller, nearing the order of magnitude of 
the effective range. For example, the virtual n- 9 Li energy E nc = —1000 keV, along with the 
effective range r nc = 1.4 fm, corresponds to a scattering length a nc = —4.2 fm. 

The study by Marques et al. [50j extracted experimental values for ■>/ (r^ n ) for Borromean 
halo nuclei using the technique of intensity interferometry along with the two neutron cor- 
relation function to study the dissociation at intermediate energies of two neutrons in halo 
nuclei. The result for 11 Li, ■sj (f^ n ) = (6.6 ± 1.5) fm, is in close agreement with both the 

LO, VvInJ = (6-8 =t 1-8) f m ; an d NLO, \/ ( r% n ) = (6.3 ± 0.5) fm, calculations found using 
the two-body virtual energy reported in [481 ]. E nc = —800 keV. However, due to the large 
uncertainties in both the theoretical and experimental values for n Li, there exists a large 
range of E nc values which would produce a {r^ n ) value in agreement with the experimental 
value of Marques et al. j50|. The LO result for 14 Be, y/ (r^ n ) = (4.1 ± 0.5) fm, is about a 
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FIG. 3: The various mean square radii for 20 C as a function of the n- 18 C two-body energy with error 
bands from the theoretical uncertainty. The LO results represented by the solid black lines, with 
error bands represented by dashed lines (as calculated in Ref. 16f] ) . The NLO results represented by 
the red dotted-dashed lines with solid error bands. As input, the ti-ti scattering length d nrh — 

-18.7 

fm, the three-body binding energy = 3506 keV, and for NLO results the effective range r$ = 1.4 
fm were used. 



Fermi below the the experimental value from Ref. 50], v (^nn) ex = (5-4 ± 1.0) fm, but the 

two values are consistent within one error bar. The NLO result, (r 2 nn ) = (3.86 ± 0.09) 
fm, is even smaller and only consistent with the experimental value within two sigma errors. 
Recall, however, that the NLO results are found using only an estimate for the effective 
ranges. Our results indicate that the effective range in the n- 12 Be-interaction is below this 
estimate. Furthermore, the large uncertainty in these experimental values is indicative of 
the need for more precise measurements of the mean square distances in 2n halo nuclei. 
Also, the recent work by Orr [5l| discusses the care which must be taken in interpreting 
the experimental results of j50| . Specifically, the technique is sensitive to the population of 
states in the continuum by the dissociation process rather than being a true ground state 
measurement. For more details on this issue, see Ref. [51] . 

Turning now to the so called Samba halo nuclei 12 Be, 18 C, and 20 C, for which the n-c 
subsystem is bound, it can be seen that a positive effective range always shifts the mean 
square radii to larger values. We focus on the case of 20 C and look at the mean square 
radii over a range of E nc values. The results, using the central value for the three-body 
binding energy as input, B { ° ] = 3506 keV, were calculated to LO in Ref. 116]. We have now 



calculated these values to NLO, using the inverse of the pion mass as an estimate of the 
effective range, and the results compared to the LO case can be seen in Fig. [3j with error 
bands estimated from the theoretical uncertainty, as described above. The solid black lines 
represent the LO results, with error bands represented by the dashed lines. The dotted- 
dashed red lines represent the NLO results, with solid colored error bands. We see that the 
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FIG. 4: The various mean square radii for the Efimov excited state of 20 C as a function of the n- 18 C 
two-body energy with error bands from the theoretical uncertainty. The LO results represented by 
the solid black lines, with error bands represented by dashed lines (as calculated in Ref. 16)]). The 
NLO results represented by the red dotted-dashed lines with solid error bands. As input, the n-n 
scattering length a nn = —18.7 fm, the three-body binding energy = 3506 keV, and for NLO 
results the effective range ro = 1.4 fm were used. 



mean square radii to NLO show the same behavior with changing E nc as the LO case, with 
an increase in E nc leading to an increase in size. Also we see more clearly the general shift 
toward larger radii due to the effective range corrections. 

As was shown in the previous subsection, there exists one Efimov excited state in 20 C for 
E nc < 165 keV. In Fig. HJ the mean square radii for this excited state to NLO are plotted 
over a range of E nc values, together with the LO results from Ref. [l6j]. Again we see that 
the general behavior of the mean square radii as a function of the two-body energy remains 
the same when calculated to NLO. At the endpoints, signifying the critical E nc where the 
excited state breaks up, the mean square radii diverge as the excited state is destroyed and 
the particles fly apart. For small values of E nc , there is a positive shift in the mean square 
radii for positive effective range, although this shift is relatively small, as compared to the 
relative shift in the ground state radii. However, as E nc grows, there is a region where a 
positive effective range causes a negative shift in the mean square radii. This is related to 
the region where the effective range causes a positive shift in the binding energy, as described 
in the previous subsection. 

In Table HJ we have highlighted the results using two different values of the n-c two-b ody 



binding energy: E nc = (530 ± 130) keV [45| and the older value E nc = (162 ± 112) keV [39 
Moreover, we show results for E nc = 60.0 keV which would lead to an Efimov excited state. 
The rows marked by 20 C* give the results for the excited state. The three-body binding 
energy of this excited state is listed in the second column with theoretical uncertainty as 
calculated in Sec. IIII Al 
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Within the error estimates our results generally agree with the values obtained previously 
by Yamashita et al. using the renormalized zero-range model 14j . Note, however, that this 
model does not include the effective range corrections and no error estimates are given. 

Finally, we note that we explicitly do not apply our method to 6 He. The n- 4 He-interaction 
is dominated by strong P-wave resonance in the J = 3/2 channel, while our effective the- 
ory assumes that S-wave interactions dominate. A treatment of 6 He therefore requires an 
effective theory with a different counting scheme. While various schemes to treat such P- 
wave resonances in Effective Field Theory have been developed [l2, T^, their application to 
three-body systems remains to be worked out. 



IV. CONCLUSION 

In this paper, we have studied the next-to-leading order (NLO) effective range correc- 
tions to the universal properties of two-neutron halo nuclei. Our investigation was performed 
within an effective quantum mechanics framework using cutoffs of the order 500 MeV. As- 
suming that the halo nuclei have resonant S-wave interactions between the neutron and the 
core, the effective potential reduces to a separable S-wave potential. At NLO the potential is 
determined by the neutron-neutron and neutron-core energies and effective ranges plus one 
three-body parameter that can be fixed from any three-body datum. At next-to-next-to- 
leading order in the expansion in Mi ow / 'M^gh no new two-body parameters enter. Whether 
another three-body parameter enters at this order is still an open question that has not been 



completely resolved [19|, [21 



Our main focus was in investigating NLO corrections to the results for In halo nuclei 
which arise from a non-zero effective range. Because the neutron-core effective range is 
not known experimentally, we were only able to estimate the range corrections assuming the 
effective range to be of the order of the inverse pion mass, r ~ l/fn n = 1.4 fm. Based on this 
assumption, we found that the leading order results are very robust and range corrections 
are typically small. We investigated the possibility of excited Efimov states in halo nuclei 
assuming that the ground state is also an Efimov state. In particular, we found that the 
corrections to the parametric boundary curves, within which at least one excited Efimov 
state will occur were essentially unchanged for realistic values of the effective range. We 
found that the value of the critical n- 18 C two-body energy required for such an excited state 
to exist in 20 C was decreased slightly from 165 keV at LO to 161 keV at NLO. While the 



TUNL nuclear data evaluation from 1995 [38l. l39| quotes this energy as E nc = (162 ± 112) 
keV, the most recent evaluation AME2003 [3 finds E nc = (580 ± 100) keV which excludes 
an excited Efimov state in 20 C. We also calculated the NLO shift in the excited state binding 
energy, finding that |A£?3^| < 0.5 keV, a nearly negligible shift from the LO result. 

Moreover, we calculated the shift from the LO results in the mean square radii for known 
2n halo nuclei. We found that for Borromean halo nuclei, the general tendency is for a 
positive effective range to shift all mean square radii to smaller values, unless the two-body 
n-c virtual energy is very close to threshold. The opposite was found to be true for the 
so called Samba halo nuclei, for which the n-c subsystem is bound. The positive effective 
range shifts the mean square radii of these halo nuclei to larger values. We also compared 
our results for the n-n mean square radius to the experimental data for the Borromean halo 
nuclei n Li and 14 Be 5(j. While there is good agreement in the n Li system, the results for 



14 Be only agree within two sigma errors. This indicates that the effective range in the n- 12 Be- 
interaction is smaller than our naive estimate tq = 1.4 fm. For a better description within our 
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effective theory, it will be necessary to determine the as yet unknown parameters, specifically 
the various neutron-core scattering lengths and effective ranges. Such information will also 
provide a solid basis for evaluating their universal behavior and determining the relevance 
of the Efimov effect as a possible binding mechanism for halo nuclei. 

If the Efimov effect is responsible for the binding of at least some three-body halo nuclei 
or their excited states it would open the possibility of a universal binding mechanism for 
higher-body halo nuclei as well. In Ref. [52|, a pair of universal tetramer states associated 
with every Efimov trimer was predicted for the case of identical bosons and their binding 
energies were calculated in the vicinity of the unitary limit. Recently, the spectrum of 
these states was mapped out for all scattering lengths and their experimental signature in 



4-body recombination was pointed out [53|. Soon after this prediction, the loss resonances 
from both tetramers in an ultracold gas of 133 Cs atoms were indeed observed [54|. Another 
calculation indicated that even higher-body cluster states with universal properties might 



exist [55[ , leading to an intriguing paradigm for a universal binding mechanism of weakly 
bound nuclei near the drip lines. Effective theories represent an ideal tool to explore this 
scenario in nuclei and predict its signatures in the structure and reactions of halo nuclei. 
Together with new precise experimental data from facilities such as FAIR and FRIB this 
will open the possibility to test this scenario. We note that this scenario does not require 
that the ground state of a halo nucleus is an Efimov state. More general situations where 
only excited states have universal character are also possible. 
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Appendix A: Wigner Bound 

In this appendix, we derive the Wigner bound on the S-wave effective range for our 
effective potential. In the renormalization process, the bare coupling constants must be 
redefined twice in order to absorb high-energy effects proportional to positive powers of the 
cutoff parameter, which is equivalent to adding counter terms to the effective potential. The 
bare couplings {C ,C 2 } of Eq. (TT]) are related to the redefined couplings {C ,C 2 } by: 



C 2 A 5 



°o = ri^TW + z-ir2n 2 2fiCi (Al) 

and 



c*= n ~„ 2Aa 1-1 + J , mr*inj \> (A2) 



■2- 



where \i is the two-body reduced mass. One clearly sees that there are real values of Ci for 
which the bare couplings become complex. This would imply a complex effective potential. 
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Therefore, we must place a constraint on the value of C2: 



= 2 A 3 



7T 3 



< 1. 



(A3) 



The renormalization of the effective potential is completed by tuning the redefined cou- 
pling constants to reproduce the scattering length a and the effective range r according to 
the following expressions: 



1 9 
= 2vr 2 2u 

C 



and 



Co 



-/(x ) - -A 
a 7r 



f /(So) - fx 



(A4) 



4 - 22 HQ/(-o)-|A)- 



(A5) 



where the parameter xo is analogous to x n i from Eq. f JTOj) found with a n ; = a and r ni = r . 
Substituting Eq. (IA5j) into Eq. (IA3j) . we find the Wigner bound for our effective potential: 



2 A 3 f/(xo)-H 
" 3 Q/(*o)-fA) : 



< 1. 



(A6) 



If the cutoff is much larger than the momentum scale of the T-matrix pole, A ^> y/2fiE 2 , 
then f(x ) — > 1 and we have: 



,2 4 

n> < 2 ~ T 



6 



7T 3 

71 A aA 2 + 2^A7 



(A7) 



It is then obvious that for A — > 00, the effective range must be equal to or less than zero. 
Conversely, by setting the input values of the scattering length and effective range, Eq. (1A6j) 
places a constraint on the maximum value of the cutoff A. 



Appendix B: Wave function and matter density form factors 

The full bound state wave function of a two-neutron halo nucleus can be reconstructed 
from the solutions for the spectator functions F n and F c found from the coupled integral 
equations (TJIE]). The form of the wave function, i{pq\^) = ^i(p, q), depends on the choice 
of two-body subsystem and corresponding spectator particle, where the index i = n,c labels 
the spectator particle. Recall that in the wave functions, the Jacobi momentum p describes 
the relative momentum between the two particles in the chosen two-body subsystem, while 
q describes the momentum of the spectator particle relative to the center of mass of the 
two-body subsystem. 

To reconstruct the S-wave part of the full wave functions with either a neutron or core 
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spectator from the spectator functions we find: 



*n(p,g) = [G n (p } q;B 3 



H(A) 
A 2 



t n {q; B 3 )F n {q) 



+ ^ dx t n (n' nn ,B 3 )F n (n' nn ) + t c (n' nc ; B 3 )F c (n, 
H(A) 



ncl 



(Bl) 



A 2 



t c (q;B 3 )F c (q) 



dx t n (n' cn ,B 3 )F n (jr' cn ) 



-i 



(B2) 



where the two-body T-matrices t n and t c are given in Eqs. (0 [8]). The expressions for the 
propagators Gq and Gq are: 



B 3 ) 

G c ( P ,q;B 3 ) 
Finally, the shifted momenta are: 



A + 1 2 A + 2 ., 



-i -l 



2A * 2(A + 1] 



q 



B 3 +p 2 + 



A + 2 
AA ' 



-i -i 



(B3) 
(B4) 



<.n = n' nn (p,q) = Jp 2 + 



2pqx, 



(B5) 




A 2 A 

< c = ^Lfe 9) = \\v 2 + ( A+l y I 2 ~ ~A+1 2P<1X ' ^ 



and 



tt',, = K l Jl>-<l)= \!p 2 + ^q 2 -pqx. (B7) 

The three-body wave functions can be used to calculate other low-energy properties of 
the three-body bound state. With the Jacobi momentum states it is straightforward to 
calculate the Fourier transform of the one- and two-body matter densities with respect to 
the momentum transfer squared. These are defined as the one- and two-body matter density 
form factors ^(k 2 ) and J r n i(k 2 ), respectively, where i — n,c. The derivation of the form 
factors from the three-body S-wave wave functions can be found in Refs. 0, Q. The 
resulting expression for the one-body form factors is: 

^{k 2 ) = Jdpp 2 Jdqq 2 J dx ^(p, q)%(p, V V + k 2 - 2qkx), (B8) 

while the two-body form factors are: 

Tnc{k 2 ) = Jdpp 2 Jdqq 2 J dx * n (p, q)y n (Vp 2 + k 2 - 2pkx, q), (B9) 
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and 



T nn {k 2 ) = Jdpp 2 jdqq 2 J dx * c (p, q)V c ( ^p 2 + k 2 - 2pkx, q). (BIO) 
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